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Abstract 

in 

0^ We discuss some theoretical uncertainties in the calculation of the cross section for 

charm production in charged current deep inelastic neutrino scattering related to ambigu- 
ities in the treatment of terms which are singular in the limit of a vanishing charm mass. 
In particular we compare the so-called variable flavour scheme where these terms are 
absorbed in the parton distribution functions containing the charm as an active flavour, 
with the so-called fixed flavour scheme with no charm mass subtraction where the charm 
appears only in the final state of fixed-order scattering matrix elements. Using available 
parametrizations of parton distribution functions we find that the two schemes lead to 
largely differing results for separate structure functions whereas the differences cancel to a 
large extent in the total cross section in that kinematical region which has been measured 
so far. 
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1 Introduction 



Information on parton distributions comes mostly from totally inclusive lepton scattering on 
nucleons. The increased precision in the experimental measurements of deep inelastic scattering 
(DIS) and related processes has lead to a considerable improvement in our knowledge of the 
parton distributions of the proton and the neutron. Various global analyses of the data have 
been done in recent years to obtain increasingly better parton densities [1, 2, 3]. All these anal- 
yses, however, give only limited information on the densities of heavy quarks, strange, charm 
and bottom, since they make a much smaller contribution to the measured structure functions. 

This situation has improved recently for the strange-sea density. To extract this density 
from deep inelastic data two approaches have been used. In the first method one considers the 
difference of the F 2 structure functions measured in neutrino and muon DIS which selects the 
strange-sea contribution [4]. This method suffers from experimental normalization uncertain- 
ties since data from different experiments are combined. In this paper we are interested in the 
second method which consists in measuring the charm production in charged current neutrino 
and antineutrino DIS where the characteristic signature is the production of dimuons in the 
final state [5]. The dimuon data obtained by the CCFR collaboration were originally analyzed 
in the simple parton model as scattering of strange quarks in the nucleon supplemented with 
the slow-rescaling prescription to account for the effect of the finite mass of the charm quark 
[6]. The result of this analysis in leading order QCD (LO) was a strange-sea contribution with 
a strength reduced by a factor of ~ 0.4 as compared to the non-strange quark and antiquark 
components of the sea and a fitted charm quark mass m = 1.31 GeV . Later this analysis 
was repeated using the next-to-leading order (NLO) formalism of Aivazis et al. [7] with the 
result that the strange-sea distribution was still reduced by a factor of 0.5 compared to the 
non-strange-sea distribution, but the charm quark mass was now changed to m = 1.70 GeV . 
This 50 % suppression in relation to the u and d sea distributions is now incorporated in all re- 
cent parton distribution fits [1, 2, 3]. Thus the s distributions in MRS(A) [1] and CTEQ3M [2] 
agree quite well with the NLO strange-sea distributions of the CCFR collaboration [5] whereas 
in the recent GRV analysis [3] the generated s distribution matches the LO distribution in the 
earlier CCFR analysis [6]. 

It is well known that in NLO the parton distribution are ambiguous and have no direct 
physical meaning since they depend on the subtraction or factorization scheme. Thus, in order 
that the s distribution assumed in the MRS or the most recent CTEQ analyses are consis- 
tent with the CCFR results, the same scheme has to be chosen. The CCFR analysis is based 
on the NLO formalism of Aivazis et al. [7]. Their calculation of the charm production cross 
section is performed in the MS scheme with subtraction of the collinear singularity m — > 0. 
This is the variable flavour scheme in the terminology of ref. [7]. This scheme is preferred for 
situations with Q 2 ^> m 2 . The factorization scale in [7] is chosen to be // = 2pj ax where 
p-max _ A(W 2 , m 2 , M 2 )/2W is the maximum available transverse momentum of the final state 
charm quark for given kinematic variables x and Q 2 . M is the nucleon mass and W the total 
cm. energy. This scale // is somewhat larger than Q depending on the relevant value of x which 
is (x v is) = 0.15 in the CCFR analysis. Except for this special choice of the factorization scale, 
which is also used for the scale of the QCD coupling constant, the variable flavour scheme cor- 
responds to the subtraction applied in the MRS and CTEQ analyses. This prescription implies 
that the charm contribution is introduced as an active flavour of the nucleon which absorbs 
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the charm mass collinear singularity in the NLO current-gluon fusion contribution. This also 
implies that the leading logarithms due to the charm mass singularity are resummed through 
the usual QCD evolution of the charm distribution. Actually, in the MRS and CTEQ analyses 
the charm mass is equal to zero except for the thresholds at Q = m introduced in a s (Q 2 ) 
and in the charm parton distribution. The GRV collaboration, on the other hand, in a recent 
structure function analysis [3], have taken the point of view that the charm contribution in the 
nucleon can be generated dynamically through perturbative contributions and with no extra 
non-perturbative charm contribution that would have to be determined from DIS and other 
experimental data. This is the so-called fixed flavour scheme with Nj = 3 in the terminology 
of [7]. Here the structure functions F% and F[ in the neutral current process are provided by 
fixed-order perturbation theory from the LO fusion process *y*g — > cc. The scale // in a s (jj 2 ) 
and the gluon distribution function is chosen as // = 2m, irrespective of Q 2 } which according to 
[8] leads to results similar to the corresponding NLO results of Laenen et al. [9]. A comparison 
of the fixed and variable flavour schemes has been performed for neutral current heavy quark 
production in [10]. 

Actually in this GRV analysis [3] the strange sea distribution is also generated dynamically, 
and the strange sea at their input scale //q, which is very small, vanishes. It is clear that in 
the fixed flavour scheme, where the charm quark is totally extrinsic, the charm production 
in charged current processes, when considered in NLO, must be calculated also in the same 
scheme. Such a calculation has been performed some time ago by two of us [11] with the inten- 
tion to study the charm mass dependence in relation to the higher order QCD corrections. At 
that time no attempt was made to analyze the scheme dependence and to compare the fixed 
flavour scheme with Nj = 3 with the variable flavour scheme with Nj = 4. Since now these 
two schemes have been applied in the construction of parton distributions for the nucleon, it 
is appropriate to find out, whether in both of these two schemes the dimuon production cross 
sections in charged current neutrino and antineutrino DIS can be obtained in NLO in agree- 
ment with experiment. Since the measurements are performed in an intermediate Q 2 range 
(~ (Qli s ) = 25.5 GeV 2 [5]) there is a priori no reason to prefer the variable flavour scheme as 
was done in the analysis of the experimental data [5]. Thus it would be advisable to repeat the 
experimental analysis with the fixed flavour scheme to obtain the NLO strange sea distribution 
in a consistent way in this scheme, too. Due to the experimental details involved in such an 
analysis this could be done only by the CCFR collaboration. Instead we take the following 
approach. The cross section for charm production depends in general on three independent 
structure functions F^ c } F^ and F^ c . We calculated these structure functions in the variable 
flavour scheme with the same input as in the CCFR analysis and consider the result as a repre- 
sentation of the experimental results, of course only in the range 0.01 < x < 0.2, as measured 
in the CCFR experiment. For the proton (neutron) structure functions we use the CTEQ3M 
or the MRS(A) parton distribution functions which agree with the CCFR NLO strange quark 
distribution function up to small deviations. This can be done easily by converting the for- 
malism in [11], which was obtained for unsubtracted charm, to the variable flavour scheme, i.e. 
with intrinsic charm. Then we calculated the F( c , F^ and F^ c structure functions in the fixed 
flavour scheme using the parton distribution functions of GRV [3] , which are the only available 
parametrizations of this scheme, and compared the results in the two schemes. In addition we 
have investigated the effect of scale variations, in particular by replacing the scale by Q 2 } which 
might be the more natural choice instead of (pj ax ) 2 . 
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The outline of this work is as follows. In section 2 we repeat some necessary formalism and 
specify our input. The numerical results for the charged current structure functions F( c , F^ 
and F^ c in the two schemes are presented and compared in section 3. Here we also discuss to 
what extent differences for the structure functions obtained for the two schemes might have 
shown up in the CCFR data, since in the CCFR experiment only some particular average over 
the three structure functions could be measured. 



2 Formalism and Input 

The signature for the production of charm quarks in neutrino- and antineutrino-nucleon scat- 
tering is the presence of two oppositely charged muons in the final state. For neutrino scattering 
on protons the underlying process is a neutrino interacting with a s or a d quark, producing a 
charm quark that fragments into a charmed hadron. The charmed hadron's semileptonic decay 
produces a second muon of a charge opposite to that of the first one. The analogous reaction 
with an ingoing antineutrino proceeds through an interaction with a s or d antiquark, again 
leading to oppositely charged muons in the final state. 

The inclusive cross section for charm production in v^N collisions (N stands for proton or 
neutron) v^N — > fi~cX is described by the structure functions F( c , F^ and _F 3 SC , which depend 
on x and Q 2 } and is given by the following formula: 



= fr 5 ( (1 " y)Fr{x ' Q2) + y 2xF "( x > Q2 ^ + y{l ~ yi 2 ) xF ^ ^ 2 )) • C 1 ) 

For the reaction v^N — > fi + cX the structure functions F-j 50 , _F 2 SC and F^ c are replaced by 
F^,TT and (-Pf ). G F is the Fermi coupling and x = Q 2 /2Pq, S = (k + P) 2 , y = Pq/Pk = 
Q 2 1 xS with P and k being the momenta of the ingoing nucleon and neutrino, respectively. The 
effect of the W propagator is neglected as usual. 

At leading order (LO), charm is produced by scattering directly from s and d quarks in the 
proton (s and u quarks in the neutron). For scattering off an isoscalar target the cross section 
is given by (I) where the structure functions in LO are given by: 

FZ c (x, Q 2 ) = x (\V cd \ 2 (u(x } Q 2 ) + d(x, Q 2 )) + \V CS \ 2 2s(x, Q 2 )) = 2xs'(x, Q 2 ) 

Fr(x,Q 2 ) = \Ft(x,Q 2 ) = s'(x,Q 2 ) (2) 

where x is the slow rescaling variable, which takes into account the non- vanishing charm mass 
m, x = (I + m 2 IQ 2 )x [12]. V c d and V cs are the CKM mixing matrix elements and s'(x) is the 
'effective' strange sea distribution in LO. The LO expression, given by (1) and (2), shows the 
sensitivity of the process to the strange quark sea distribution. Charm (anticharm) production 
from scattering off d (d) quarks, respectively u (u) quarks, is Cabibbo suppressed. In the case 
of charm produced by neutrinos, half of the cross section is due to scattering from s quarks, 
although the d quark content of the proton is ten times larger. In antineutrino scattering, 
where d (u) quarks from the sea contribute, approximately 90 % of the cross section is due to 
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scattering of s quarks. In next-to-leading order, we have contributions proportional to 0(a s ) 
which come from virtual corrections to the LO diagrams, gluon emission contributions and con- 
tributions from VK-gluon fusion. The LO and NLO hard scattering diagrams are shown in Fig. 
1. Due to the size of the gluon distribution, which is an order of magnitude larger than the sea 
quark distributions, the gluon-initiated diagram is of similar magnitude as the LO contribution. 
The NLO quark-initiated diagrams in which a gluon is radiated are smaller, since they are not 
enhanced by a large parton distribution. 

In our earlier work [11] (see also [13]) we have calculated the 0(a s ) corrections to the struc- 
ture functions Ff c } F^ and F^ c for a finite c quark mass m on the basis of the diagrams in Fig. 
1 using dimensional regularization. The formulas given in [11] are for the fixed flavour scheme 
with the DIS convention for the subtraction of the collinear singularity of the intermediate s 
quark. In this scheme the F^ c structure function, for example, was given by (see (2.3b) in [11]) 



" 2 ) = s'(x,Q 2 ) + 



ol s{i j2 r) f 1 dz 



2tt 



+ 



2 \ 

fl9 (*) 



d Z — 

u qg z, n2 



9[-,Q'' 



(3) 



In the following we shall define the various subtraction schemes for the contribution induced 
by the gluon in the nucleon, i.e. by the process W + g — > cs which is the term proportional to 
g(x/z,Q 2 ) in (3). 

In the scheme with no charm subtraction the functions d qg (z, ^y) and f qg (z) have the fol- 
lowing form: 



m 



d qg = T R \(1- z) 2 + z 2 + iz— (1 - z 



n9 9 m 2 / m 2 
:i-zY + z 2 + ^—[l-2z-z— 



hv 



^2 \ 



+ [(l-zY + z 2 } In 



(1 ~ z, 

z(l-z) 



(4) 




+ 8z(l-z)-lj (5) 

where s = Q 2 (l — z)/z and Tr = 1/2. In (4) we must distinguish z and z = (1 + m 2 /Q 2 )z. In 
[11] the factorization scale for the s quark was chosen as fj 2 = Q 2 . If we have an arbitrary scale, 
it appears as an additional term in (5). The expression for f qg (z) in (5) corresponds to the DIS 
scheme. In (5) we corrected an error [14] which is already present in the original calculation 
[15]. The f qg given in [11, 15] is not exactly in the MS scheme because initial gluon spins were 
averaged in four rather than in n dimensions. The form as given in (5) must be used for the case 
that the neutrino cross sections are calculated with nucleon structure functions constructed in 
the DIS scheme. When using structure functions in the MS scheme, (5) is replaced by the 
simpler expression 

^(^) = ^((l-^) 2 + ^ 2 )ln^. (6) 



5 



In the variable flavour scheme, the leading singularity in the charm mass oc In m 2 is sub- 
tracted as well. For the MS scheme, this subtraction is very simple. We combine it with d qg} 
which then contains all mass dependent terms. Then the total d^i s is 



= d qg + A<- (7) 

with 



A^ = r,((l-,-) 2 + ,- 2 )ln^ 



where fi s is the subtraction scale for the charm term. We see that the term Ad qg subtracts just 
the contribution proportional to lnm 2 in (4), and the total d^i s has no singularity for m — > 0. 



In the DIS scheme we must subtract additional terms which are equal to the contribution 
to the F 2 structure function in neutral current deep inelastic scattering originating from the 
transition ^g — > cc. This subtraction is not unique. In our calculation we keep the charm 
mass m finite whereas in all structure function analyses (except those of GRV) the c mass 
is neglected except for a threshold in the Q 2 evolution related to m. From this viewpoint it 
would be reasonable to subtract terms corresponding to the massless cross section for ^g — > cc 
except for the term proportional to ln(/i 2 /m 2 ). In this case the subtraction term has the same 
structure as f qg in (5) and d qg has to be modified by an additional term AJ^ 15 which is 

Ad° IS = T R (j(l - zf + z 2 ] (hi + In ^pj + 8*(1 - z) - l) . (9) 

Of course, (8) still contains mass dependent terms since z is replaced by z, i.e. we use the 
slow-rescaling form of the variable z as in f qg . This is partially dictated by the requirement 
that the terms proportional to lnm 2 must cancel [13]. We observe that with this prescription 
terms proportional to m 2 In m 2 survive. 

Another possibility is to follow Gottschalk [13] and to subtract the contribution of F 2 for 
massive quarks which is obtained from ~fg — > cc with m ^ 0. Here we have the problem that 
the thresholds for ~fg — > cc and W + g — > cs are different, 4m 2 for the first process and m 2 for the 
second process. The subtraction in [13], which we call the DIS m scheme, proceeds as follows. 

2 2 

The cross section for jg — > cc is given by c 9t 2(z, ^y) (this is the term equivalent to d qg (z, ^y) 
in (3)) [16]: 



R 



z) 2 + z 2 + z(l 



3z) 



4m 2 



,8m 4 



In 



1 + /3 



8^(1 - z) - 1 - z(l - z) 



1-/3 
4m 2 



(10) 



In (10), z is the variable with no slow-rescaling factor, i.e. z = Q 2 /{Q 2 + s) and 



0=1 



4m z 

Q 2 l -z 



1/2 



4m 



2\ 1/2 



11] 
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The slow-rescaling variable lor ^g — s- cc is ~z = (Q 2 + 4m 2 )/(Q 2 + s) = (1 + 4m 2 /Q 2 )^, which 
has the property, that lor J — > 1 we have s — > 4m 2 , so that the argument ol the square root 
in (11) is always positive. This is not the case il we would express z by J and would consider 
z — > 1. The subtraction term in the DIS m scheme can be written as: 

AC 5 " = (*, ^) § + Mil - =zf + t] In i£ (12) 

where z must be expressed by the variable ~z using the relation given above. In a second step 
~z is identified with the slow-rescaling variable z Irom the process W + g — > as. This procedure 
is certainly not unique. (12) differs Irom (8) by terms proportional to m 2 /Q 2 . For Q 2 ^> m 2 
these two ways ol subtraction give approximately the same result. 

We shall not write down the equivalent subtraction terms lor the contributions induced 
by the s quark in (3), i.e. by the process W + s — > eg. In the scheme with no charm subtrac- 
tion the lunction d qq (z } ^y) was given in [11] together with the subtraction term f qq (z) in the 
DIS scheme lor fj 2 = Q 2 . The form ol this term for fj 2 ^ Q 2 is obvious Irom (5) and the 
MS subtraction term is analogous to (6). The charm subtraction terms AJ^ 5 , AJ^ 15 and 
Ad q J q ISm are constructed analogously to (7), (9) and (12). In the same way one proceeds for the 
structure lunctions F^ c (x } Q 2 ) and _F 2 sc (x, Q 2 ). For these the 0(a s ) corrections have been given 
also in [11] and the various subtraction terms are constructed in the same way as for F^ c (x } Q 2 ). 

The calculation ol the structure lunctions F-^ c (x, Q 2 ), F^ix, Q 2 )-, F^ c (x, Q 2 ) and 

F[ c (x } Q 2 ) = FZ c (x, Q 2 ) - 2x F lS c(x, Q 2 ) (13) 

for charm production by neutrinos in the various schemes has been done with the following 
nucleon structure lunction sets. For the variable flavour scheme, where the charm is treated 
as intrinsic, we have chosen the sets CTEQ3 [2] and MRS(A) [1] both in the MS and DIS 
schemes. For these we have chosen the scale = fi s = 2pj ax and m = 1.7 GeV as in the 
CCFR analysis [5]. In the following we shall refer to this scheme also as the Nj = 4 scheme. 
In the fixed flavour scheme, also denoted as the Nj = 3 scheme, we work with the parton dis- 
tributions GRV(94) [3]. In this new structure lunction analysis the charm quark distribution 
to the electromagnetic structure lunctions F\ and F 2 is generated dynamically Irom *y*g — > cc 
at the scale = fi s = 2m with m = 1.5 GeV . These parton distributions exist also in a MS 
and a DIS version. To obtain the structure lunctions F t sc ( they appear in the CCFR 

analysis we have calculated them for an isoscalar target including the Cabibbo suppressed terms 
ol the charged current with \V c d\ = 0.224. This value is to be compared with the PDG value 
V us = 0.2205 ± 0.0018 and with \V cd \ = 0.232 + 0.018(-0.020) as determined in the CCFR 
analysis [5], lurthermore |Ks| 2 = 1 — |Kd| 2 - 

In the next section we present our numerical results and compare the charm production 
structure lunctions for CTEQ, MRS(A) and GRV with each other for the MS and DIS schemes. 

3 Numerical Results 

Before we come to the charm production structure lunctions, we compare the most important 
input, i.e. the s quark parton distribution ol the proton. It is plotted at fj 2 = 25 GeV 2 for 
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CTEQ3M and MRS (A) in the MS version and GRV(94) also in the MS version as a function 
of x between x = 0.01 and 1 in Fig. 2 including the Cabibbo suppressed contribution coming 
from u and d quarks for the case of an isoscalar target, which is (see (2)): 

s'(x, ii 2 ) = \V cs \ 2 s(x, ii 2 ) + \V cd \ 2 -(u(x, ii 2 ) + d(x, ii 2 )). (14) 

We see that the CTEQ3 and MRS (A) distributions are almost equal in the considered range 
of x whereas the s' distribution of GRV(94) differs from the other two. It is clear that this 
difference comes from s(x). The small admixture from u(x) and d(x) in (13) has no effect since 
the valence parts of these quark distributions are equal to those of MRS (A) by construction [3]. 
Since GRV use a different scheme for the charm quark, there need not be any agreement for 
the s quark either since the distribution functions in NLO have no physical meaning. The only 
requirement is that the usual deep inelastic and other data are reproduced. We also mention 
that the strange quark distribution xs(x) in MRS(A) [1] has been compared to the distribution 
obtained in the CCFR analysis [5]. At fj 2 = 4 GeV 2 it lies inside the error band for xs(x) as 
given by CCFR and so satisfies the experimental constraint of reproducing the dimuon data if 
calculated with the same NLO formalism. 



Next we present the comparison of results for the structure functions in the MS scheme. 
We choose Q 2 = 25 GeV 2 since this corresponds very well to the average (Q 2 ) = 25.5 GeV 2 in 
the CCFR analysis. Figs. 3, 4, 5 and 6 show the plots for xF{°(x, Q 2 ) } F£ c (x } Q 2 ) } xF^ c (x } Q 2 ) 
and F[ c ( function of x in the interval [0.01, 1]. We observe that for all four structure 

functions there is little difference between the CTEQ3 and MRS (A) results as to be expected 
from the comparison in Fig. 2 (note that also the gluon densities of the CTEQ3 and MRS (A) 
parametrizations are very similar). If we compare with the LO approximation which is es- 
sentially given in Fig. 2, up to a factor of 2 for F^ and xF^ c } we see that for CTEQ3 and 
MRS(A) the 0(a s ) corrections together with the subtraction terms change xF^ c } F^ and xF^ c } 
in particular towards small x. xF^ c and F^ have decreased and xF^ c has increased strongly. 
In the GRV scheme, xF( c and F^ increased as compared to xs'(x), but xF^ c now decreased 
strongly. One can understand this pattern by looking at the different subtraction terms for the 
two cases which were given in the last section. The dominant subtraction term for the Nj = 4 
structure functions (CTEQ3 and MRS(A)) contains the factor ln(/i 2 Jm 2 ) with /i 2 s = i(pj ax ) 2 , 
p-max _ (yy2 _ m 2^ and rises like ln(l/x) for decreasing x, whereas the equivalent factor for 
the GRV distribution, ln(Q 2 /m 2 ), is independent of x. Thus, the increase of the GRV strange 
sea towards small x is reflected unchanged in the results for the structure functions, whereas 
the less steep strange sea of the CTEQ3 and MRS(A) distributions is made even less steep. 
We emphasize that the scales in the Nj = 4 and in the Nj = 3 case are different. For x > 0.1 
the functions xF^ c and F^ are almost equal in the two schemes but they differ appreciably for 
x < 0.1, whereas xF^ c is much smaller in the Nj = 3 scheme for all x. F[ c shown in Fig. 6 
differs also in the two schemes, but now more for x > 0.05. However F[ c is small and certainly 
very difficult to be measured. 

In case that the three structure functions could be determined experimentally by measuring 
the charm production cross section in (1) as a function of y with fixed x and Q 2 , the behaviour 
of these structure functions as a function of x, as shown in Figs. 3, 4 and 5, could be verified. It 
is clear that the structure functions in the two schemes are not compatible with each other, so 
that either one of them or both must be adjusted. This could be obtained by a modification of 
the input s'; but one would have to check also whether a change of the input gluon distribution 
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would lead to a better agreement without coming into conflict with other measurements. Un- 
fortunately the cross section (f ) could be measured only in a rather limited range of y values. 
In the CCFR analysis (y) ~ 0.47. Therefore it makes sense to compare an effective structure 
function F^ for this fixed value of y. As in [11] we introduce 

F^{x, Q\ y) = (l- y)F?{x, Q 2 ) + y 2 xF[ c (x, Q 2 ) + y 2 xF* c (x, Q 2 ), (15) 

which is just the expression in the curly brackets in (1), and evaluate it for y = 0.5 in the two 
schemes. The result is plotted in Fig. 7. Comparing the results we see that the F^ are not 
so different. The differences are below 20%. Below x = 0.03 F^ in the GRV scheme is larger 
by up to 15% and for x > 0.03 it is smaller by up to 20%. Since xF( c and F^ on one side 
and xF^ c on the other side influence F^ in opposite directions these differences cancel each 
other to a large extent (see Figs. 3, 4, 5). Whereas, for example at x = 0.01, the ratios r of the 
results for the GRV and CTEQ3M parametrizations for F^ and xF^ c are r = 2.3 and r = 0.36, 
respectively, the ratio is decreased to r = 1.15 in F^ ■ One also notes that the difference in 
xF^ c is responsible for the reduction of F^ in the GRV scheme for x > 0.1. These differences 
are significant even when we consider that the measurements of CCFR have an average experi- 
mental error of 10% as one can see from their results for the s quark distribution. Therefore, we 
might expect that, if the CCFR analysis would be repeated in the fixed flavour scheme with the 
GRV structure functions as input, the prediction of the charm production cross section would 
deviate from the measured data. By how much is difficult to assess without going through the 
analysis, which must include all the experimental details like cuts, for example. 

In the fixed flavour scheme the subtraction terms are quite different from those in the vari- 
able flavour scheme, i.e. we have very different 0(a s ) corrections. In addition, a much smaller 
factorization scale was chosen in the former. Therefore also xs(x) has to be different in the two 
schemes. From Fig. 7 it appears that xs(x) in the GRV set can be adjusted in such a way that 
the differences in the 0(a s ) correction terms for F^ between the two schemes can be compen- 
sated. Whether this can be done for the separate structure functions is an open problem. The 
figures show that xF( c and F^ are still approximately equal to the structure functions in the 
variable flavour scheme in the region x > 0.1, but differ appreciably for x below 0.01, whereas 
xF^ c differs considerably in the whole x range. As an exercise we repeated the calculation 
in the fixed flavour scheme with xs'(x) taken from MRS(A) and xg(x) taken from GRV(94). 
The differences for F^ between the two schemes were reduced for large x but remained at 
low x and there was no improved agreement for the separate structure functions. A sensible 
test to distinguish the two schemes would be the separation of the three structure functions by 
measuring the y dependence of the charm production cross section. This y dependence reflects 
itself in F^ in the following way. For example at y = 0.3, the ratio r is increased to 1.4 at 
x = 0.01 or, at y = 0.7 and x = 0.1, we find r = 0.76. 

We studied the structure functions in the variable flavour scheme also with the scale Q 2 
instead of (2pj ax ) 2 . With this scale xFf c (x, Q 2 ), F^ix, Q 2 ) and xF^ c (x, Q 2 ) look somewhat 
different, but F^^x^Q 2 ) changes very little over the whole range of x > 0.01. Thus consid- 
ering this effective structure function as the only quantity measured so far, the CTEQ3M and 
MRS(A) parton distributions account very well for the CCFR data also with the scale Q instead 
of 2pj ax . Results for the scale dependence of individual structure functions are presented in 
Figs. 8, 9, 10, 11 and 12. In these plots we show the structure functions for the MRS(A) parton 
distributions using the scales = fi s = f-2p^ ax with / = 1/2 and / = 2, respectively (dashed 
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lines, results for the CTEQ parametrization are similar). The scale changes are significant only 
for small x. Here xF^ c shows a somewhat stronger variation than the other structure functions 
_F 2 SC and xF^ c . In F^ at y = 0.5 (Fig. 12) this somewhat stronger variation is reduced again, 
so that the predictions are rather reliable with respect to the CCFR cross section. This agrees 
with similar findings in [5]. To obtain an estimate of the scale variation also for the GRV case 
we replaced the scale 4m 2 by AQ 2 . The results, contained also in Figs. 8-12 (dotted lines), 
show the same pattern: individual structure functions are slightly decreased (increased for F^ c ), 
however F^ shows only small changes. 

We also made the comparison in the DIS scheme. The results are essentially the same. The 
resulting structure functions in the GRV scheme are almost equal to those in the MS scheme 
since the parton distributions are changed only by the appropriate subtraction terms which 
are compensated by identical subtraction terms for the charm structure functions. The same 
applies for the MRS (A) parametrizations whereas for CTEQ3D new fits have been made and 
therefore the differences are slightly larger. It turns out, however, that the calculated charm 
structure functions xF( c , F^ and F^ c are again almost the same. In the DIS scheme the differ- 
ence in F^ between the Nf = 3 and the Nf = 4 schemes is even smaller than it was for the 
MS scheme. However, the differences in xF( c , F^ and F^ c between the GRV and the variable 
flavour scheme stay essentially the same also in the DIS scheme. We checked that the two 
possible treatments of the subtraction in the DIS scheme as given in Eqs. (9) and (12) lead to 
small shifts of the results for xF( c , F^ and F^ c reaching a few percent only at small x, but the 
influence on F^ is negligible. 



4 Concluding Remarks 

In this article we have discussed the scheme dependence of predictions for the cross section 
of charm production in charged current deep inelastic neutrino scattering. We compared two 
different schemes, the so-called variable flavour scheme with 4 active flavours where the charm 
mass singularity is absorbed in the parton distribution functions, and the fixed flavour scheme 
with only 3 active flavours and no charm mass subtraction. This comparison could be performed 
in a consistent way by using corresponding parametrizations of parton distribution functions 
for N f = 4 and for N f = 3. 

We have seen that the two schemes lead to largely differing results for separate structure 
functions, but the differences cancel to a large extent in the effective structure function F^ in 
the kinematical region of intermediate values of y ~ 0.5 which has been measured in the CCFR 
experiment so far. 

It is a general consensus that the scheme dependencies are not physical. This means it 
should be possible to adjust in either scheme the input distribution functions, in particular 
xs(x) and xg(x) } in such a way that physical observables for charm production are reproduced. 
Taking the point of view that the structure functions with Nf = 4 and CTEQ3 or MRS (A) 
as input represent the "experimental data", it would be an interesting exercise to find strange 
quark and gluon distributions which would reproduce these "experimental data" in the Nf = 3 
scheme. Of course, it would be preferable to start from real experimental data for the separate 
structure functions which unfortunately do not exist yet. 
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ure Captions 



. 1 Diagrams contributing to charm production in neutrino scattering: a) Leading order 
diagram and gluon initiated diagrams, b) Diagrams with gluon radiation and virtual 
corrections. 

. 2 Cabibbo-rotated strange quark distribution xs'(x,/j 2 ) (Fl c,LO (x, /i 2 ) in leading or- 
der) for /i 2 = 25 GeV 2 in the ITS parametrizations of CTEQ3M (full line), MRS (A) 
(dashed line) and GRV(94) (dotted line). 

. 3 Structure function xF( c (x) for charm production in neutrino scattering in the MS 
scheme using the parametrizations of CTEQ3M (full line), MRS (A) (dashed line) 
and GRV(94) (dotted line). 

. 4 Same as Fig. 3 for the structure function _F 2 sc (x). 

. 5 Same as Fig. 3 for the structure function xF^ c (x). 

. 6 Same as Fig. 3 for the structure function F[ c (x). 

. 7 Same as Fig. 3 for the structure function F^^ (x) (see eq. (15)) for y = 0.5. 

. 8 Scale dependence of the structure function xF( c (x) for charm production in neutrino 
scattering in the MS scheme using the parametrizations of GRV(94) with scale = 
2m (upper dotted line) and = 2Q (lower dotted line) and for the parametrizations 
of MRS(A) (dashed lines) with scale = fi s = f ■ 2pj ax . Upper line: f = 1/2, 
lower line: f = 2. 

. 9 Same as Fig. 8 for the structure function _F 2 sc (x). Upper lines: = 2m and / = 1/2, 
lower lines: = 2Q and / = 2. 

. 10 Same as Fig. 8 for the structure function xF^ c (x). Upper lines: = 2Q and / = 2, 
lower lines: = 2m and / = 1/2. 

. 11 Same as Fig. 8 for the structure function F[ c (x). Upper lines: = 2m (low x) and 
/ = 2, lower lines: = 2Q (at low x) and / = 1/2. 

. 12 Same as Fig. 8 for the structure function F^^ (x) for y = 0.5. Upper lines: = 2m 
and / = 2, lower lines: = 2Q and / = 1/2. 
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